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The "optical springs" regime of the signal-recycled configuration of laser interferometric 
gravitational-wave detectors is analyzed taking in account optical losses in the interferometer arm 
cavities. This regime allows to obtain sensitivity better than the Standard Quantum Limits both 
for a free test mass and for a conventional harmonic oscillator. The optical losses restrict the gain 
in sensitivity and achievable signal-to-noise ratio. Nevertheless, for parameters values planned for 
the Advanced LIGO gravitational-wave detector, this restriction is insignificant. 
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I. INTRODUCTION 

Signal-recycled "optical-springs" topology of interfero- 
metric gravitational- wave detectors [ll I2| is now consid- 
ered as a likely candidate to design the second generation 
of these detectors, such as Advanced LIGO Ij. This 
topology offers an elegant way to overcome the Standard 
Quantum Limit (SQL) for a free mass — a characteris- 
tic sensitivity level when the measurement noise of the 
position meter is equal to its back-action noise 0, 0, 0] ■ 

This method is based on the use of optical ponder- 
motive rigidity which exists in detuned electromagnetic 
cavities ^ ^] . It turns the test masses in a gravitational- 
wave detector into harmonic oscillators, thus providing a 
resonance gain in test masses displacement signal [lOllllj . 
In the signal recycling topology, the pondermotive rigid- 
ity can be created relatively easy by adjusting position 
of the signal recycling mirror. In this case, only the 
signal (anti-symmetric) optical mode eigenfrequency is 
changed. The arm cavities and the power (symmetric) 
mode remain resonance-tuned causing, therefore, no ad- 
ditional problems with the pumping power. 

In large-scale optical systems having bandwidth equal 
to or smaller than the signal frequency f2, the op- 
tical rigidity has a complicated frequency dependence 
0, IT^ [3j ■ This feature allows to obtain not one but 
two mechanical resonances and consequently two minima 
in the noise. Alternatively, these minima can be placed 
close to each other or even superimpose, thus providing 
a single wider "well" in the noise spectral density. 

The former reg ime was exp lored in detail in several 
articles P, 0, Q, E, E, 0|- The latter was exam- 
ined rather briefly in papers 0, 0|. At the same 
time, it looks very promising for detection narrow-band 
gravitational-wave signals with known frequencies, in 
particular, those from the neutron stars. In this article 
we present an in depth analysis of this regime. 

In Sec.^ we analyze dynamic properties of the 
frequency-dependent pondermotive rigidity in no optical 
loss case, including, in particular, an instability inherent 
to the electromagnetic rigidity. This section is based, in 
part, on the results obtained in the articles [l5| 

In Sec. lIIII we compare the optical rigidity-based 
scheme with other methods of circumvent the SQL for 
a free mass and show that it has a significant advantage, 
namely, it is much less vulnerable to optical losses. 

In Sec. lIVI we calculate the sensitivity of the Ad- 



vanced LIGO gravitational- wave detectors in the double- 
resonance optical springs regime. 

To clarify this consideration we try to avoid bulky cal- 
culations in the main text. In Appendix 1X1 we give de- 
tailed analysis and calculations of sensitivity and signal- 
to-noise ratio for Advanced LIGO interferometer. In Ap- 
pendices |b1 and El we provide calculations for simphfied 
model without optical losses. 



II. FREQUENCY-DEPENDENT OPTICAL 
RIGIDITY. NO OPTICAL LOSSES 



"Conventional" v.s. "double" resonances 



In a single Fabry-Perot cavity, the pondermotive rigid- 
ity is equal to flUT^ : 



v{n) ^ [-iVL + + 5'' , (1) 



where $7 is the observation (side-band) frequency, ujp is 
the pumping frequency, 5 = ojp — ujo is the detuning, 
LOo is the cavity eigenfrequency, L is its length, £ is the 
optical energy stored in the cavity, and 7 is the cavity 
half-bandwidth. 

In articles 0, 0| the signal-recycled gravitational- 
wave detectors topology was considered in detail and its 
equivalence to a single cavity was shown. It was shown, in 
particular, that Eq. ^ is valid for this topology too, with 
obvious substitution of 7 and 5 by the anti-symmetric op- 
tical mode half-bandwidth 70 and detuning (5o [see Ap- 
pendix |^ and Eqs. HA26p . (|A28|l ]. Energy £ in this case 
is equal to the total optical energy stored in both inter- 
ferometer arms: 



£ = 



(2) 



where Ic is the power circulating in each arm. For the 
sake of consistency with other parts of the paper, nota- 
tions So and 70 will be used throughout this paper. 

It should be noted that for the narrow-band regimes 
which we dwell on, 70 have to be small: 



70 < (^0 . 
Therefore, we neglect for a while term 70, setting 
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(3) 



(4) 
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where Kq is the rigidity value at zero frequency: 



(5) 



A detailed analysis is provided in III Bl and IIVI 

Consider a harmonic oscillator having mass m and 
rigidity (Q. We can write the equation of motion as 
follows: 



m—^ + m6l'^^^^+5lKox{t) 



dPF{t) 



5lF{t), 
(6) 

where x is the oscillator coordinate and F is an external 
force acting on it. 

In general the system has two resonances with frequen- 
cies 



^4 



± 



(7) 



In articles P, 0, they are referred to as "mechanical" 
and "optical", because in asymptotic case Kq 0, the 
frequencies are equal to: 



(8) 



Therefore, the high-frequency resonance can be readily 
interpreted as a result of the optical power sloshing be- 
tween the optical cavity and detuned pumping field, and 
the low-frequency one — as a conventional resonance of 
mechanical oscillator with Kn rigidity. 

The second-order pole [Ij, or double resonance 
case takes place if these two frequencies are equal to each 
other: 



v2 



i.e. if 



where 



Kn 



ttiSq 



4^ 



crit : 



£rrit 



8ujQ 



(9) 



(10) 



(11) 



is the critical energy. In this case, the equation of motion 
© has the following form: 



d^Fjt) 
dt^ 



2nlF{t) 



(12) 



It is useful to consider action of resonance force F{t) 
Fq cos rtot on this system. Solving Eq. 112|l , we obtain: 



xit) 



Ik 
8m 



-t cos Hot 



t sin ripi 



(13) 



The leading term in amplitude of x{t) grows with time 
as . At the same time, for a conventional oscillator we 
have 



x{t) 



FqI sin Qnt 



of- 



(14) 



and for a free mass, 



x{t) 



Foil 



cos 



Hot) 



(15) 



Therefore, response of the "double resonance" oscillator 
on the resonance force is (ftot/A) times stronger than that 
of a conventional harmonic oscillator, and (fio^/^)^ times 
stronger than that of a free mass one. 

It was shown in 0,^3' that due to this feature, the 
"double resonance" oscillator has much smaller value of 
the Standard Quantum Limit for narrow-band signals, 
than both free mass and conventional harmonic oscilla- 
tors. It is convenient to express this gain in terms of 
dimensionless parameter 



(16) 



where Sh is the single-sided spectral density of detector 
noise normalized as the equivalent fluctuational gravita- 
tion wave h, 



8h 



(17) 



is the value of Sh corresponding to SQL. For conventional 
SQL-limited gravitation wave detectors (with free test 
masses), C ^ 1- 

It is shown in Appendix ^ that in case of a conven- 
tional first-order resonance, the equivalent noise curve 
has a "well" at resonance frequency flo, which provides 
the gain in sensitivity: 



_ An 

?osciU — \l 



(18) 



where AD, is the bandwidth where this gain is provided. 
In particular, this gain can be obtained by using either 
the "mechanical" or "optical" resonance of optical rigid- 
ity)- 

In the case of "double resonance" oscillator, this gain 
can be substantially more significant: 



An An 
^'^"-T^ ^^"^^"VTi^- 



lis) 



Even better result can be obtained if optical energy is 
slightly smaller than the critical value (|ll|l : 

Ko^^il-i]^) ^ f = £:crit(l-r7'), ?y«i- 

(20) 

In this case, function ^(fi) has two minima, which corre- 
spond to two resonance frequencies 



r!± = r!oVT±^«f^o(i±| 



(21) 



and a local maximum at frequency f2o- H parameter rj is 
equal to the optimal value 



(22) 



then the bandwidth is ^/2 times wider than in pure dou- 
ble resonance case for the same value of ^ : 



dbl 



?cnhdbl = 



An 



V2 V2rjo 



(23) 
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Above considerations give an important result concern- 
in^ the signal-to-noise ratio values for different regimes 

The signal-to-noise ratio is equal to: 



SNR: 



Shin) 



dn. 



(24) 



where h(Q) is the gravitation wave signal spectrum. It is 
shown in Appendix [CI that for a conventional resonance- 
tuned interferometer (without optical springs), 



01\-[\no springs A 



'SQL 



(25) 



where factor Af for wide band signal (with bandwidth 
Arisignai > ^^o) is about unity: for short pulse Af w 2.0 
and for step-like signal A/" « 0.7. For narrow-band signals 
with Afigignai <C we have Af ~ Ar^signai/f^o- 

In the narrow-band cases described above the signal- 
to-noise ratio can be estimated as follows: 



SNR' 



\h{no)\^An 



An 



|fe(»o)P 



(26) 



(we omit here a numeric factor of the order of unity). 
In case of a conventional harmonic oscillator, it can be 
shown using Eqs. I|18I26|I . that: 



SNRosciii 



2|fe(»o)P^^o 



(27) 



(numeric factors in this formula and in Eq. H28|l below 
are obtained through rigorous integration in Eq. (|24|) . 
see Appendix [nj. This value does not depend on Afl. 
Thus conventional harmonic oscillator provide arbitrary 
(in the case of zero losses) high sensitivity at resonance 
frequency, and also gain in signal-to-noise ratio equal to 
~ ilo / Ailsignai for narrow-band signals as compared with 
conventional interferometer. At the same time, the value 
(|27|l is fixed and can not be increased by improving the 
meter parameters, and there is no gain in signal-to-noise 
ratio for wide-band signals. 

On the other hand, in case of double-resonance oscil- 
lator, it stems form Eqs. (|19l I25|l . that: 



SNR, 



dbl 



V2a 
An 



\h{no)\^no 

X 



no 



^sql(^o) 



V2An 



X SNR„ 



(28) 

Therefore, the double resonance oscillator allows to in- 
crease both the resonance frequency sensitivity and the 
signal-to-noise ratio by decreasing the bandwidth An. 

The signal-to-noise ratio for Advanced LIGO interfer- 
ometer with optical rigidity and influence of optical losses 
is considered in Sec. lIV 



B. Dynamic instability 

It is well known 's','^ that positive pondermotive rigid- 
ity is accompanied by negative dumping, i.e. a ponder- 
motive rigidity-based oscillator is always unstable. In 
case of frequency-de pen ded rigidity, this instability was 
calculated in article (see also 0,0). 



If the eigenfrequencies are separated well apart from 
each other, n+ — ri_ ^ 70, then the characteristic insta- 
bility time equals to 



^instab 



70^0 

ni-ni 



7o 



0.1-=- Is. 



(29) 



The instability becomes more strong, if n^ — n_ 0. 
In double resonance case (51+ = fi-), the instability time 
is equal to 



^instab 



Viono 



(30) 



However, inequality f^o'Hnstab S> 1 holds in this case too. 

It has to be noted that, in principle, any instability 
can be dumped without affecting signal-to-noise ratio, if 
the feedback system sensor sensitivity is only limited by 
quantum noises but its implementation is a separate 
problem which we do not discuss here.. 



III. COMPARISON WITH OTHER METHODS 
TO OVERCOME SQL AND INFLUENCE OF 
OPTICAL LOSSES 

A. "Real" vs. "virtual" rigidities 

From the Quantum Measurements Theory point of 
view, a laser interferometric gravitational-wave detector 
can be considered as a meter which continuously moni- 
tors position x{t) of a test mass m 0,0,0 • Output signal 
of this meter is equal to: 



i{t) ^x{t) + £niotor(t) , 



(31) 



where imotor(i) is the measurement noise and x{t) is the 
"real" position of the test mass. It includes its responses 
on the signal force 



-^signal (t) 



mLh{t) 



and on the meter back-action force -F'motcr(0- 

X{t) = Z-l[Fsignal(t) + ^^mctcr(i)] 



(32) 



(33) 



(the terms containing the initial conditions can be omit- 
ted, see article 7]). Here Z is a differential operator 
describing evolution of the test object. For a free test 
mass {i.e. for the initial LIGO topology). 



Z = m 



dt^ 



(34) 



If rigidity K (including the frequency-dependent ponder- 
motive one) is associated with the test mass, then 



(35) 



Therefore, the signal force estimate is equal to 

F{t) = Zi{t) = Fsignal(i) + Fsum noise (t) , (36) 
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where 

-^suin noise (^) — -^mctcr(^) "t' ^•^mctcr(^) ■ {^^) 

is the meter total noise. 

The back-action noise i^inctcr(0 proportional to the 
amplitude quadrature component of the output light. 
The measurement noise Xmctcr(0; the simplest case 
of SQL-limited detector, is proportional to the phase 
quadrature component. In this case the spectral densities 
of these noises satisfy the following uncertainty relation 
(see e.g. Q): 

SxSf > fi^ ■ (38) 

In more sophisticated schemes which allow to overcome 
SQL, the noises -Fmctci(i) and Xmctci(i) correlate with 
each other. In this case, the back-action noise can be 
presented as follows: 

-^inctcr(0 — -^mctcr(0 ~^ ^*^mctcr it), (39) 

where ^m°ctcr(0 the back-action noise component non- 
correlated with the measurement noise and coefficient JC 
can be referred to as "virtual" rigidity. 

Note that it is precisely the idea of quantum variational 
measurement 0, |23, ll^, ll^l • In conventional opti- 
cal position meters, including the LIGO interferometer, 
one measures the phase quadrature component in output 
wave. This component contains both the measurement 
noise (a;mctcr) produced by phase fluctuations in input 
light wave and the back action noise (-Fmctcr) caused by 
amplitude fluctuations (optimization of the sum of these 
two uncorrelated noises produces SQL). However, using 
homodyne detector one can measure tuned mix of the 
phase and amplitude quadratures of output waves and 
this mix can be selected in such a way that the back 
action noise can be compensated by the noise of the am- 
plitude quadrature. It can be considered as introduction 
of correlation between the back action and the measure- 
ment noises as presented in Eq. (|39|l . We see that the 
"virtual" rigidity only relates to the measurement proce- 
dure (homodyne angle). 

Substituting Eq ^ into Eq. (jSH, we obtain, that 

-^sum noise (0 — -^meter(0 ^~ ^cff'-^mctcr (0 i (^^) 

where 

Zeff = Z + /C = m— + K,s , (41) 

and 

K,s = K + IC. (42) 

is the effective rigidity. 

Thus, in the lossless case, the total meter noise 
only contains the sum of real rigidity K and virtual one 
IC, and replacement of any one of them by another one 
does not change the total noise spectral density and the 
signal-to-noise ratio |24j . 

Spectral densities of the noises i;mctcr(i) and F^^^^^.{t) 
also satisfy the uncertainty relation 

^.4°^ > , (43) 



which does not permit simultaneously making both noise 
terms in Eq. (|4()|l arbitrary small. However, factor Zcff 
can be made equal to zero by setting 

JsTcff = mfl^ . (44) 

In this case, only noise F^^^^j.{t) remains in Eq. (|40() . In 
principle this noise can alone be made arbitrary small, 
thus providing arbitrary high sensitivity. 

Both "real" and "virtual" (created by the noise correla- 
tion) can be used for this purposes. Applying the idea of 
variational measurement a simple frequency-independent 
cross-correlation (and thus frequency-independent vir- 
tual rigidity /C) can be created relatively easily by using 
a homodyne detector. In this case, Eq. (|44|l is fulfilled at 
some given frequency, creating resonance gain in sensitiv- 
ity similar to one provided by a conventional harmonic 
oscillator. 

A frequency-dependent cross-correlation and thus a 
frequency-dependent virtual rigidity K. can be induced 
through modification of the input and/or output optics 
of the gravitation-wave detectors using additional large- 
scale filter cavities In this case, condition (|44|l can 
be fulfilled and thus the sensitivity better than SQL is 
obtained in any given frequency range. 

Consider now the real pondermotive rigidity. It is ev- 
ident that it can not be tuned in such a fiexible way as 
the virtual one. In the double-resonance case, graphics 
of K{^1) and mil^ touch each other only at frequency Hq. 
However, around this point, 

{K,s{n)) - ^{n- flof , (45) 

instead of 17 — fio for the case of a conventional frequency- 
independent rigidity. As it was shown in Sec.m slightly 
better results can be obtained by using sub-critical pump- 
ing £ < ScTit , which provide two closely placed first-order 
resonances, see Fig.^b). 

B. Influence of optical losses 

It is possible to conclude from above considerations 
that the virtual rigidity provides a more promising so- 
lution than the real pondermotive one. However, as 
we show in this subsection, the pondermotive rigidity 
has one important advantage: it is much less vulnera- 
ble to losses in optical elements. Our consideration will 
be based on the following statement, which follows from 
Eas. (|A39|I : a lossy optical position meter is equivalent 
to the similar lossless one with gray filter attached to its 
signal port. This filter transmittance has to be equal to: 

_,load 

■ = — — (46) 

-^equiv ; \^^J 

70 

where 

7o=7r'+7r, (47) 

the term describing signal mode half-bandwidth 
7o describing coupling with photodetector and 7q°'^'' is the 
term describing optical losses. 
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FIG. 1: Plots of the frequency-dependent rigidity KiyC) as a 
function of V? (solid line); dashed line corresponds to mfi^; 
(a) — critical pumping (the double resonance regime); (b) — 
sub critical pumping. 



Term "similar" means, that (i) both meters have same 
dynamic parameters, i.e. bandwidths and eigenfrequen- 
cies of the cavities, etc. In particular, half-bandwidth of 
the lossless meter has to be equal to 70 and (ii) optical 
energies stored in the lossless meter and the lossy one 
have to be equal to each other. 

Therefore, the lossy meter output signal can be pre- 
sented in the following form (compare with Eq. (|31|l ): 



ilossy(0 = X(t) + imetor(0 + [A l^loss (i) 



where 



l-4o| 



i-r2 . 

cquiv 



cquiv 



_ .load 

7o 



(48) 



(49) 



is the effective loss factor and iioss(i) is the additional 
noise produced by the optical losses uncorrelated with 
ammeter- Spectral density of this noise is also equal to S^^ 
because both iioss (0 and iioss (i) are in essence zero-point 
fluctuations normalized by the transfer function of mea- 
surement system. The key point is that while the back- 
action noise i^mctcr(i) can correlate with the measurement 
noise i:motor(i), it can not correlate with additional noise 

■^loss (^) ■ 



(t) + \Ao\ZS:ioss(t) 
(50) 

[compare with Eq. (pHl ]. 

It follows from this equation that in the lossy meter 
case there is no symmetry of the real and virtual rigidi- 
ties: the term proportional to the meter noise a;motor(i) 
still depends on the effective rigidity, but the new losses 
term only depends on the real rigidity. 



sum noise V 



Compare now too particular cases of "pure real" and 
"pure virtual" rigidities. 
In the flrst case, Zes = Z, and 

-^sum noise 

(t) + \Ao\ ■^loss it)]- 

(51) 

Therefore, both terms proportional to Xmeter(0 and 
S^ioss(t) can be canceled by setting Z = 0, thus providing 
arbitrary high sensitivity at least for one given frequency. 
In the second case, 

^sumnoise(i) = ter (^) + Zeffimeter (i) ""^^^^ |A l^loss (0 ■ 

(52) 

Suppose that Zoff is canceled using some QND procedure. 
In this case, the sum noise will still consist of two non- 
correlated parts proportional to F^ctcrit) and xiossIOj 
and its spectral density will be equal to: 



(53) 



Taking into account uncertainty relation (|43(l . it is easy 
to see that 



/ >5'sum(^) 



ssQ^m 



loss 

7g 

_,load 

7o 



1/4 



where 



(54) 



(55) 



(56) 



[see Eqs. (HZEll)]. 

The restriction (|55|1 shows that the use of "pure vir- 
tual" rigidity (variational measurement) is very sensi- 
tive to losses — for parameters from Table we have 



(^loss/^load 



1/4 



0.7 (!). 

For conventional (resonance-tuned) Advanced LIGO 
topology, 70 ~ ~ 27r X lOOs"^ > 7^°'''' and "virtual" 
rigidity can be introduced by means of the variational 
measurement. In this case the optical losses restrict the 
sensitivity by the following value: 



loss 

To 



1/4 



0.2 



(57) 



(for the value of 7i°^^ see TableHJ. 

In general both real and virtual rigidities exist in 
the signal-recycled configuration of the intcrferomctric 
gravitation-wave detectors as well as for that of a single 
detuned cavity. However, in the narrow-band regimes 
(PI, the real rigidity dominates: it is approximately 6q/jo 
times stronger than the virtual one (compare Eqs. (jA34|) ') 
and HA50p , or (jA53ap and (jA53b|) , which differ by terms 
of the order of magnitude < 7o/(5o only). Due to this 
reason, this regime is free from limitation H55|l . 

It is also interesting to note that while both K and 
JC contain, in general, imaginary parts, in effective rigid- 
ity Kcs these imaginary parts exactly compensate each 
other. 
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TABLE I: Parameters planned to use in Advanced LIGO Q- 
In estimates of \Ao\'^ and 70 using formulas l|4'7l4i)IA26t we 
assume that 11 + 7?«e^"*P ~ 2 



Q = 0, 77 = 0, 70/^0=0.01 



Transmissivity of SR mirror 


rri2 
^ S 


= 0.05 


Transmissivity of input mirrors in 


rp2 _ 


= 0.005 


arms 






Loss coefficient of each mirror in arms 


Al- 


= Al = 1.5 X lO"'^ 


Length of interferometer arm 


L = 


4 km 


Effective loss factor 


\M\ 


2 = 0.24 


Relaxation rate of difference mode 


70 = 


= 2.9 s"^ 


"Intrinsic" relaxation rate 


„,loss 

7o 


= 0.56 s-i 


Mean frequency of gravitational wave 


no = 


= 27r X 100 s"^ 


range 







IV. ADVANCED LIGO SENSITIVITY IN 
DOUBLE-RESONANCE REGIME 

A. The sum noise spectral density 

Sensitivity of the signal-recycled Advanced LIGO 
topology in the narrow-band approximation (close to the 
double resonance) is calculated in Appendix ^ It is 
shown, that in this case [see Eq. (|A57p ] 



loss (~) 3 



where j/ = 17 — 



0; 



loss 

7o 
^o\Ao.\ 



■C, C = l 



V2 



(58) 



(59) 



and ?7q, \Aa\ are basically parameters 77, |^o| corrected 
to take into account virtual rigidity (i.e. homodyne angle 
a), see Eqs. (|A55|) . 

It follows from Eq. H58|) that the single minimum de- 
pendence of sensitivity ^{fl) takes place if 770, = with 
bandwidth Afl equal to 



An 



'27i°-f7o 



(60) 



On Figs .[3 El we present the sensitivity plots for various 
sets of parameters which allows to realize single mini- 
mum dependence of S, on different frequencies Qq and 
with different depths for the case when effective loss fac- 
tor \Aq\'^ — 0.24 (as planned in Advanced LIGO, solid 
curves) and for no losses case (dotted curves). We see 
that the sensitivity degradation due to the optical losses 
is negligibly small. Possibility of scaling the frequency 
ilo is also demonstrated in this plots. 

For r] > we have two minima dependence of sensitiv- 
ity ^ at two different frequencies n± ~ ^q^/I ± 77. With 
increase of 77 the distance between minima also increases. 
Comparing values ^(ri±) with ^(r2o) we can introduce a 
"characteristic" value rjac when V2^{fl±) — ^{flo): 

Co 



In this case 



P - 



loss 



^0\Aa 



c An^ 



"0 



(61) 



(62) 
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FIG. 2: Top plot: function v'Sh(r2)//isQi(noo) plotted for 
parameters a = 0, 70/^^0 ~ -01, \Aaf = 0.24 and -qa = 0. 
The left solid curve — Qo = 0.5 x Qqo, the middle solid curve 
— ilo = 1^00 and the right solid curve — f2o = 2 x Qoo- The 
dotted curves — sensitivity for no loss case with 70° = 
70. Straight line depicts SQL sensitivity. Bottom plot: the 
same as on top plot with homodyne angle a = n/2. Qoo is 
some arbitrary frequency, for example, Q,oo — 2it x 100s~^. 




1 

FIG. 3: Left plot: sensitivity curves ^(fi) plotted for parame- 
ters a — 0, 0,0 = 0,0, \Ao\'^ = 0.24 and 77 = 0. The wider solid 
curve corresponds to 70/^0 ~ -01, the narrower solid curve 
— 7o/f^o ~ .001. The dotted curves — sensitivity for no loss 
case with 70° '""""^ = 70. Straight line depicts SQL sensitiv- 
ity. Right plot: the same as on the left plot for homodyne 
angle a = n/2. 



[compare with Eq. (|23|l ] . The sensitivity plots for this 
case (again for loss factor |-4o|'^ — 0.24) are given in 
Fig.|31 

On Fig. |S1 we also present sensitivity curves for well 
separated minimums (77 = 277c) and slightly less relax- 
ation rate 70 = O.OSfio fot different values of homodyne 
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a = 0, T] ^ r]c 



a = Tv/2, T) = T]c, 




[compare with Eq. H57|) ]. Here for the estimate we used 
parameters from Tabled 

Unfortunately, this value can not be obtained for 
planned Advanced LIGO parameters. Indeed, it follows 
from Eqs. (|X26l) . l|X28|) . that 



So < 



27^ 



load 



r2 



(64) 



FIG. 4: Left plot: sensitivity curves ^ plotted for parameters 



(The physical sense oi'^}^'^'^ /T^ is quite clear — it is relax- 
ation rate of a single FP cavity in one arm.) In the double 
resonance regime we need to have ~ \/2^o ~ 10^ s~^. 



0, 7o/f2o = .01, \Aof = 0.24 and = r]c- The dot- However for Advanced LIGO parameters (Tabled 



ted curves — sensitivity for no loss case with 70° '"^^^'^ = iq. 
Straight line depicts SQL sensitivity. Right plot: the same as 
on the left plot with homodyne angle a = n/2. 



a = 0, 7t/2 




1 

FIG. 5: Sensitivity curves ^ plotted for parameters 70/^^0 ~ 
0.03, \Ao\'^ = 0.24, 77 — 2r)c, a — Q (curve with wider sepa- 
rated minimums) and a — 7r/2. The dotted curves — sen- 
sitivity for no loss case with 70 ° '"'^'^'^'^ — . Straight line 
depicts SQL sensitivity. 



angle a = 0, 7r/2. It is close to regime considered in Ij 
including the negligible role of optical losses (compare 
with plots on Fig. 8 in \^). 



B. Minimum of the spectral density 

Suppose that almost monochromatic signal has to be 
detected and we are interested in the minimum of nor- 
malized spectral density ^ at some fixed frequency. First 
of all we have to remind that 70°'''' is fixed (it depends on 
mirror's absorption only) whereas the effective loss factor 
\Ao\'^ (as well as |-4aP) can be modified by variation of 
mirrors transmissivities. 

It follows from Eq. H58() that ^(fi) reaches its minimum 
at J/ = ±ria^Q/2 and \Aa\ = 1, and this minimum is 
equal to 



3 

V2 



0.06 



(63) 



HI 



load 



T2 



100 s" 



(65) 



It is less by one order of magnitude than required. 
The problem can be solved through modifying Advanced 
LIGO parameters, namely, by decreasing the signal recy- 
cling mirror transmittance by approximately one or- 
der of magnitude and corresponding increase in the arm 
cavities input mirrors transmittance by the same value. 



C. Signal-to-noise ratio 

As was mentioned above, the double-resonance regime 
allows to obtain signal-to-noise ratio better than SQL 
even for wide-band signals. It was shown in [l^ that for 
no loss case the gain in signal-to-noise ratio, in princi- 
ple, can be arbitrary high. Here we show that optical 
losses restrict this gain. The details of calculations are 
presented in Appendix IA 91 

As an example of wide band signal we consider the 
perturbation of metric having the shape of a step function 
in time domain and the Fourier transform equal to 



h{il) = const/il. 



(66) 



It is worth to underline that the result practically does 
not depend on the shape of wide band signal spectrum 
(as alternative example one could consider a short pulse 
(delta- hmction) — its Fourier transform is a constant). 

To demonstrate the gain we also calculate the signal- 
to-noise ratio SNRconv for conventional LIGO interfer- 
ometer |23| without signal recycling mirror with registra- 
tion of phase quadrature and take the quotient V of SNR 
by SNRconv in order to characterize the gain in signal- 
to-noise ratio (the value SNRconv we calculated numeri- 
cally): 



SNR 



SNRc 



SNRc 



0.7 X 



(67) 



Using the accurate formulas (|A44|I we numerically cal- 
culated plots (see Fig. for gain V as function of rj 
at fixed ratio 70/^^0 = 0.0046 and loss factor factor 
\Aq\'^ = 0.24 corresponding to Advanced LIGO parame- 
ters (see Table U). We see that the degradation due to 
losses is not large as compared with no loss case. 

The gain V decreases with increase in rj bigger than op- 
timal — i.e. when double-resonance (77 = 0) transforms 
to two well separated first-order resonances (77 ^ ?7c). 



a = 0, -yo/^Q = 0.0046 



V 




0.1 0.2 n 0.3 

a = 7r/2, 7o/!:!o = 0.0046 



0.4 



0.5 




FIG. 6: Numerically calculated plots of signal-to-noise ratio 
gain P(r?) for different homodyne angles: a = (top) and for 
a = 7r/2 (bottom). On each plot the upper solid curve corre- 
sponds to no loss case, the lower dotted curve — to planned 
in Advanced LIGO losses with factor |^oP = 0.24. All plots 
are calculated for fixed ratio 70/^^0 = 0.0046, planned in Ad- 
vanced LIGO. 



We see from plots in Fig. |H1 that there is an optimal 
value of parameter -q when gain has maximum. Analysis 
presented in Appendix I A 91 shows that gain V reaches its 
maximum for optimal values of rja and \A\a and it is 
equal to 



0.6 X 



/ »0 



20 



(68) 



where estimates given for parameters listed in Tabled 

Note that this gain may be achieved at Advanced 
LIGO parameters — in bottom plot on Fig.|H|the maxi- 
mum of dotted curve is quite close to T^max- 

It is worth noting that in the "pure" double-resonance 
regime (77^ = 0) the gain in the signal-to-noise ratio only 
slightly differs from the maximum gain: 



' max ^ 



1.14 



(69) 



Using the signal-to-noise ratio for a conventional oscil- 
lator 1)2 7|l . one can calculate the gain 



SNR„ 



SNRc 



2.8 



(70) 



Comparing H68I) and (|70|1 wc sec that "double resonance" 
regime provides gain in signal-to-noise ratio about 7 
times larger than conventional oscillator. 



No losses: |ro| = 0, r; = 0, 7^ ^"'"'""'/no = 0.01 



X 




0.8 o.y 



Optical losses: |r-o| = 1, ~ 0, 70-1-/^0 = 0.01 



X 




FIG. 7: The plots of squeezing factor X as function on fre- 
quency fl for zero losses, homodyne angle varies from a = 
(lower curve) through a — it/6, n/3 to a = n/2 (upper curve) 
and 77 = 0. Top: no loss case, 7?° '™'''7no = 0.01. Bottom: 
non-zero losses case, |ro| — 1, 70+7^0 = 0.01. 



D. Squeezing in output wave 

It is also useful to find how optical losses affect output 
wave squeezing. Using Ea. (|A43|l one can calculate the 
squeezing factor X (for the coherent quantum state, X = 
1): 



x{n) 



ikbt 



('^ vac '3' vac H~ '^-vacO^vac) 



(71) 



The result of this calculation yields squeezing factor X 
as a function of SI. For particular case 77 = the profile 
is given in Fig.d The top plot correspond to the lossless 
case (^0 = 0); the bottom one — to the case when \ Ao\ = 
1. 

We see that in both cases the squeezing monotonously 
increases with increase in homodyne angle. More im- 
portant is the fact that the values of squeezing factor are 
close to one. It confirms our assumption that it is the op- 
tical rigidity rather than pondermotive nonlinearity, i.e. 
the meter noises cross-correlation (as source of squeez- 
ing) that produces the major input into the sensitivity 
gain. 
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V. CONCLUSION 

The optical rigidity which can be created in the 
signal-recycled configuration of laser interferometric 
gravitational-wave detectors turns the detector test 
masses into oscillators and thus allows to obtain narrow- 
band sensitivity better than the Standard Quantum 
Limit for a free test mass. This method of circumnebting 
the Standard Quantum Limit does not rely on squeezed 
quantum states of the optical filed and due to this it is 
much less vulnerable to optical losses. 

Moreover, sophisticated frequency dependence of this 
rigidity makes it possible to implement the "double res- 
onance" regime which provides narrow-band sensitivity 
better than the Standard Quantum Limits for both a 
free test mass and an conventional harmonic oscillator. 

The "double resonance" regime may be useful to de- 
tect narrow band gravitational waves, e.g. from pulsars. 
Knowing pulsar parameters one can tune the bandwidth 
and sensitivity in the optimal way. It is important that 
this tuning may be produced "on line" by varying sig- 
nal recycling mirror position and adjusting circulating 
power. 

Another advantage of the "double resonance" regime 
is its better sensitivity to wide-band signals. While an 
conventional harmonic oscillator provides approximately 
the same value of the the signal-to-noise ratio as a free 
test mass, in the case of a "double resonance" oscillator 
this parameter is limited only by the optical and me- 
chanical losses and other noise sources of non-quantum 
origin. Estimates based on the Advanced LIGO param- 
eters values shown that the "double resonance" regime 
can provide more than tenfold increase of the wide-band 
signal-to-noise ratio. 
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APPENDIX A: ANALYSIS OF ADVANCED LIGO 
INTERFEROMETER. 

1. Notations and approximation 

We consider Advanced LIGO interferometer with sig- 
nal recycling (SR) mirror and power recycling (PR) mir- 
ror as shown in Fig|Sl SR, PR mirrors and beam split- 
ter are immobile and have no optical losses. We assume 
that both Fabry-Perot (FP) cavities in the east and north 
arms are identical: each input mirror has transmittance 
T < 1 and reflectivity i?i = ^/l-T^ - Aj, each end 

mirror has reflectivity i?2 = \/ 1 — , where Ai, A2 are 
optical loss coefficients for input and end mirrors cor- 
respondingly. The end mirrors and input mirrors have 
equal masses m and they can move as free masses. 




FIG. 8: Scheme of Advanced LIGO 



The contribution of losses in SR mirror and in beam 
splitter is small as compared with contribution of losses 
in mirrors in arms which is larger by factor ~ — 
consequently we assume in our analysis that SR mirror 
and beam splitter have no optical losses. 

We consider "dark port" regime when no regular op- 
tical power comes to detector in the signal port. In this 
case the power recycling mirror is used only to increase 
the mean power delivered to the beam splitter and no 
fluctuational fields from the laser (west arm) access de- 
tector in the south arm. 

Upper letters denote the mean complex amplitudes, 
lower case letters denote operator describing fluctuations 
and signal, south arm, bn is the operator of Fig|SJ). The 
electric field E in propagating wave can be written as 
sum of large mean field and small component (see details 
in 23]): 



27r fkJo 



A 



a e 



Sc \ 2TrJ 

I = hujp\A\ , a = a{uJo + cl- = cl{^o ~ ^^) 

[a(wo±r2), a+(wo±17')] = 27r (5o(r2 - 11'), 

where A is the complex amplitude, S is the cross section 
of the light beam, c is the velocity of light, / is the mean 
power of the light beam, a and are the annihilation 
and creation operators. We consider sidebands ± 
about carrier with side band frequencies f2 in grav- 
itational wave range (ri/27r G 10 . . . 1000 Hz). Detailed 
notations of wave amplitudes and mirror displacements 
are given on Fig|S| 

2. Arm cavities 

First we consider the fields propagating in FP cavity 
in the east arm shown in dashed box on Fig|Sl (formulas 
for FP cavity in the north arm are the same with obvious 
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substitutions in subscripts E N). Below we assume 
that the following conditions are fulfilled: 



< 1, T, Ai, A2 < 1, 



(Al) 



where L is the distance between the mirrors in arms 
(4 km for LIGO). 
We start with a set of equations for mean amplitudes 

BEi=iTA- RiAei, B = iTAEi~ RiAe, (A2) 
Ae\ — DEI e ' — —ti2DEi e ' . 

Below we assume that the carrier frequency of the inci- 
dent wave ojo is equal to the eigenfrequency of each FP 
cavity, i.e. e'"°^/'^ = 1. So we obtain 

Bei=TAe, Be2 — —Bei, Be — TZAe, (A3) 
iT 2i7ioad 



r 



1 - R1R2 T7 



1 - R1R2 7 ' 

7 = 7load + 7I0SS, 7- = 7load — 7loss, 
_ _ cA^ 

7load "4^' '/loss ' 



(A4) 
(A5) 



A = JAj + Al 



To calculate fluctuational components of field we start 
from the set of equations: 

bEi = iTttE ~ RiaEi + iAiCEi + RiAei 2ikyE, (A6) 
Be iTttEi - RidE - RiAe 2ikyE + iAiCei, (A7) 
aEi = OR2 {-a2E - A2E ^ikxE) + iA2e2E- (A8) 

Here eiE, e2E is vacuum fields appearing due to losses. 
We denote 9 = e*('^o+^)^/c using obvious approximation 
9 ~ l + inL/c, 

Substituting (|A8|) into (|A6() and using conditions (|A1() 
we obtain: 

, T- , T- , TaBEi2ik{xE-yE) , 
Obi ^ % flfi + e^r To H — , ( A9) 

h2E — -bEi, a2E — bEi, CLEi -bEi, (AlO) 
iT 2i7ioad 



eE 



BElAi - e2EA2 

A ' 



(All) 



Here we introduced operator ce of vacuum fluctua- 
tions, which fulfill usual commutator relation ^CEi^o + 
17), e+(wo + n')] = 2TrSo{n - fl'). 

Substitution of HA9|I into HA7|I under conditions (jAll) 
allows us to get the formula for reflected field: 

bE = a-E Ti-n - ierTTn - ATTq 2ik{xe - Ve), (A12) 
-Ri+R2{l~Al)9^ j^+in 



Tin 



1 - i?ii?26l2 



"f — ifl 



(A13) 



3. Beam splitter 

We assume that lossless beam splitter has transmit- 
tance and reflection factors equal to i/V^ and — 1/\/2, 
correspondingly. The additional phase shift 7r/2 is added 
to the west arm and — 7r/2 to the south one. It is conve- 
nient to introduce new variables 



«(±) 



Ue ± UN 



b(±) = 



bE ± bN 



V2 ' ^/2 • 

For these variables beam splitter equations read: 



(A14) 



aw — — 

a.{+) = —bw 



-1 = -b. . 



(Al5a) 
(A15b) 



These equations are valid both for the zero and the first 
approximations as the beam splitter is a linear system. 

Now we can consider the fields in entire interferometer. 
For our mean amplitudes of (±) variables we have 



B 



(±) 



UA 



(±) 



(A16) 



It is easy to note that the set of equations splits into two 
independent sets: one for W and (-I-) and second for S 
and (— ). The "S, (— )" set has only trivial zero solution 
since Ad = 0. Therefore, 

A TZA 
Ae ^ An ^ -j=, Be^ Bn ^ (Al7a) 

Aei - Ami = ^ ' , Bei =Bni = ^, (A17b) 

^From this point on we omit indices -I- for the sake of 
simplicity, A is the complex mean amplitude left in the 
beam splitter. 



4. Output field 

The signal wave, registered by detector in the south 
arm ("5, (— )" mode) is coupled with the differential 
motion of the mirrors, i.e. it has a part proportional 
to differential coordinate 



[XE - Vn) - {XN - Un) 



(A18) 



coupled with the gravitation-wave signal. We are inter- 
ested in "S", (— )" mode and below we omit subscripts (_) 
for a(_), 6(_). 

^From equations for the south arm 

bo = -RsO-D + iTs9sas , (Al9a) 
bs = -Rs9% + iTs9saD . (A19b) 

(here 6g — gi(iJo+f2)Vc^ ; jg optical length between 
SR mirror and input mirror in the arm) one can obtain 
the following formulas for ffuctuational component of the 
output field: 



bo ~TZsaD— TsTn{2ikBix + iAe), 

ce - gn 

S = 7= 

V2 



(A20) 
(A21) 
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Here we introduced the "generalized" transparency 7^ 
and reflectivity TZg of "S*, (—)" mode of interferometer: 



% = 



1 + nnRs0^s ' 



Us = 



1 + TlnRsOl 



(A22) 



Due to small value of Z <C we consider below 9s a 
constant not depending on SI. 

It is useful to rewrite formulas (|A22|) in approximation 
denoting = ie*"^: 



r 

70 



^load 



70 - iSo, r_ 



load 

7o 



loss 

-7o > 



G = 

= 70- 
70- 



1 + R^e^'-'t' 



7ioad 



load 

= 7o 
cT^T? 



loss 

7o ' 



ll + 

Tloss 



Z = 



Rse^"t'\ 

_ 

~ AL ' 
2i^^7ioad sin 2(j) 
\l + Rse^^'t>\^ 



4L|1 + i?se2*'^| 



2i?,7i°^^ sin 20 



7^2 

^ s 



(A23) 

(A24) 
(A25) 

(A26) 

(A27) 
(A28) 



Here i5o is the detuning, 70 is the relaxation rate of our 
difference mode ("S", (— )" mode) which can be presented 
as a sum of "loaded" and "intrinsic" rates. "Intrinsic" 
relaxation rate 79°''*' is exclusively provided by intrinsic 
losses in mirrors in arms while "loaded" relaxation rate 
^load pj-Qvided only by transparencies of SR and input 
mirrors. Note that for Advanced LIGO relaxation rate 
7o ~ 2 s~^ — it is much less than the mean frequency of 
gravitational wave range ~ 27r x 100 s^"^ . 



5. Pondermotive forces 

To calculate pondermotive force acting on movable 
mirrors in interferometer arms we should write down the 
equation for differential field bi = (bEi — b]\ri)/y/2 in 
approximation (|A1(I : 



[l + Rse^ 



2Bikz 
T 



(A29) 



The incident wave acts on the mirror with force pro- 
portional to square of amplitude module; and we keep 
only the cross term of this square. 



F{n) 



-int 



— , F{n) ^ hk{A*a + Aat) 



First we write down the formulas for forces acting on 
the back mirrors. The difference between the forces act- 
ing on the east and north back mirrors is equal to: 



F2 ~ 2hk {A*a2 + A2a+_ + B*&2 + ^26^ 



(A30) 



where A2 = {Ae2 - An2)/\/^-, a2 = {aE2 ~ aN2)/V2 
and so on (recall that we continue considering "S'(— )" 
mode). One can extract two terms in formula for F2: 



F2 = F„ 



-Frigid , 



(A31) 



where the first term corresponds to fluctuational compo- 
nent (back action) and the second one — to the regular 
force depending on mirrors positions (optical rigidity). 
The formula for -Fmctcr is the following 



Frr 



2ihuJoTsT BI e"f'{aD + Aoe) 
(1 + Rs e2*<^) (r - in) 

{l + Rs e-2«<^)(r* - in) 



where 



An 



A{1 + Rs e^^^ 



, (A32) 



(A33) 



is the effective loss factor. 

The equation for optical rigidity K has the following 
shape 



K = 



-Frigid 16ftA:2 7ioad'5o|-Bi| 



2? = (r-ir2)(r* -if7), u 



^ Ski Jo 
LV 



, (A34) 
(A35) 



We see that optical rigidity K is proportional to detun- 
ing Sq of the difference mode of the interferometer. This 
detuning can be only introduced by displacement of SR 
mirror (while Fabry-Perot cavities in arms remain in op- 
tical resonance). 

In approximation (|Aip the forces acting on the back 
mirrors are approximately equal (with negative sign) to 
the forces acting on input mirrors — the difference is 
negligible. So for the difference coordinate x we have the 
following equation: 



where 



and 



z{n)x = F„ 



z{n) 



lieiial : 



-mn^ + K 



F 



signal 



(A36) 



(A37) 



(A38) 



is the signal force due to action of gravitational wave, h 
— is the dimensionless gravitational-wave signal. 



6. Output signal 



We see from formula l|A32|l that the back action force 
is produced by sum of fluctuational fields: ao from sig- 
nal port and e due to losses in mirrors. It is useful to 
introduce a new pair of independent fluctuational fields 
p and q (the new basis) as following: 



p = G 



T* + in au + Aoe 



r 

AqGD 



in v/1 + IAP 



(A39a) 
(A39b) 
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These two bases are equivalent — both pairs an, e and 
p, q describe vacuum fluctuations (we do not consider 
possible squeezing of field an)- 

Then we can rewrite formulas for the fluctuational 
force Ffi and output field 6^) in a more compact form: 



SihkBl^o{l + Rse^'^) p 

i'fl ^ — X 



TTse''^{T* ^1 + I A 



TsTnBiikz. 



(A40a) 



(A40b) 



Calculating the difference position x from l)A36|l and 
substituting it into l)A40b|l one can obtain the final for- 
mula for the output field bo- 



3_D 



iJT* 



(A41) 



x{p\l- 



-P- 



—iJ-joG \ ihy/2JjoG 



fi2 (r2 + r22) ; /^^^^ -ifi) \ ' 



where 



J 



8klc 
mL 



(A42) 



This formula has two terms. The first one (~ q) is pro- 
portional to the effective loss factor and appears due 
to the optical losses. The second term has the same form 
as for no losses case — compare with Eq.(16) in nil — 



with weighted multiplier + |^op and substituted 

damping rate 70 accounting for losses in Q. It is the 
background to interpret formula (|A41|I as field reflected 
from lossless interferometer (with damping rate 70) and 
then passed through grey filter with total loss factor |^o| 
which decreases our field and adds fluctuations (the first 
term) . 

It is important that the second term has multiplier 
mn'^/Z{il) which increases the relative contribution of 
the second term at frequencies close to mechanical res- 
onance (when |^(ri)| <C mfl^) while the first term does 
not depend on frequency — it demonstrates the advan- 
tage of "real" rigidity and may explain the relatively weak 
degradation of sensitivity due to optical losses. 



7. Sensitivity 

In gravitational wave antenna one registers the quadra- 
ture component 6^ of output fields using the balanced 
homodyne scheme (not shown in FiglHl: 



V2 



(Z(f])V2(l + |AI 



{A, {qA*e-'^ + q+Aoe'"^) 



'App + A^pt+A, 



A, = 



Ar. 



A.. = 



-mn^ + ^ - 1 ■^^^ 



g-<[(r;)"+n=-iyr;-H*y7o+e='°] 

(r;)2 + n2 



no 



in T* 



in 



-y/8 J70 (70 — i^) sin a — 60 cos a 



no 

a = C-(j) 



V 



(A43a) 
(A43b) 
(A43c) 



(A43d) 
(A43e) 



Here and below we normalize dimensionless metric by 
SQL sensitivity /isQL(f^o) at some frequency ilg. 

Now we can write down one-sided spectral noise den- 
sity Sh recalculated to variation of dimensionless metric 
ha and normalize it to SQL sensitivity /isQL(f^o): 



2|A,p + 2|Ap 
\As? 



P1+P2 + P3 



(A44a) 



n^-n^iSl~j^+) + j{So-losm2a)\ , 

(A44b) 

P2 = 7o(2(5of^^- J(l-cos2a))^ (A44c) 

P3 = \Ao\'{ [n^ - (61 + -iDn" + J5oY + 47o'f^'}, 

(A44d) 



4J7or2^ 

02 

"0 



I (70 — iil) sina — (5o cos a| , (A44e) 



Note that without signal recycling mirror SQL sensitivity 
in LIGO lossless interferometer can be achieved at work- 
ing frequency fio = 7 if the optical power is equal to 
the optimal one Isql{^o) [lll23j|: 



1sql{^q) 



mn'oLc 
8lOo 



or 



J_ 



1. 



(A45) 



Although presentation (|A44|) is compact and conve- 
nient for numeric estimates, it can mask the physical 
structure of the noise. Due to this reason, we provide 
a more transparent form of this equation: 



Shin) 

where 



'-L^n^ 



n 



\v\- 



2m J Jo I (70 - in) sin a — So cos a 
is the measurement noise, 

•S'loss — \Ao\'^Sx 

is the noise created by the optical losses, 
Zcff = Kcs - mn^ , 

and 



(A46) 

(A47) 

(A48) 
(A49) 
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K,ff = 



iJ 



(A50) 



+ 70 (7o - (5q + fi^) sin 2a - 26oJo cos 2a 
is the effective rigidity. 

8. Narrow-band case 



Suppose that the observation frequency is close to the 
double resonance frequency fio = / V2 and the pump- 
ing power is close to the critical value: 



J 



V2 



where <C 1, if <C 1. 
In this approximation, 

P2~272rj[j(l + cos2 a)\ 
Q « 2^270^0(1 + cos^ a), 



(A51a) 
(A51b) 



(A52a) 

(A52b) 
(A52c) 
(A52d) 



0.2 0.4 0.6 O.S 1 1.2 1.4 1.6 1.8 2 

I -4a I 

FIG. 9: Dependence of the signal-to-noise ratio on \Ac,\ at 
optimal ria IIA6III . 



where 



C=l + -^|^„|2 + |^„|4. 



(A58) 
(A59) 



and 



9. Signal-to-noise ratio 



-7 ^ I /I 2 2r)2 ^o7o . 



V2 



sin 2a 



'5'a; 



\/2TO7of^o(l + cos^ a) 



(A53a) 
(A53b) 

(A53c) 



In this approximation we have the following formula for 
the sensitivity ^: 



70 (1 + cos^ a) 



1 



2V27of^i](l + cos2 a) 
2 



X <j (4z.2_r;2r!2-25^sin2a 



+ (4^2-»7^f^^) +4f^2^2 I (A54) 



Using the following notations: 



2 2 I 



7o sin 2a 



N/2r!o(l + |^oP) 
2_ V2\Ao\' 



1 + cos-^ a 
and taking into account that 

1 



(A55a) 
(A55b) 



70 = 7o 1 + 



(A56) 



Eq. l|X54)l can be rewritten in a more compact form: 



(A57) 



In the narrow-band case, the main contribution into 
the integral of the signal-to-noise ratio is produced in 
vicinity of ^Iq. In this case we can use approximation 
(|A57p and expand the limits of integration over from 
—00 to 00 and making substitution x = \/2iy/noy/^ C-^^^ 
with notation a ~ 'ql^/2g^/C: 



SNR: 



2 p \h{n)\'^dn 2 \h{nfi 



12 /-oo 



dv 



2\h{n)\^nQ ^ 1 



dx 



|/i(r2)|2f7o ^ 1 



1 1 

+ 



|/i(r2)|2i]o 1 VoTTPTT 



The maximum of this expression is achieved, if 
1 2 2.g%/C 

and it is equal to: 



where 



fc(l^al) 



3^/^|^o| 
2C3/4 



(A61) 



(A62) 



(A63) 



14 



is dimensionless function plotted in Fig.|5| It attains the 
maximum at 



and it is equal to 




(23 + 3\/73)3/4 



0.43. 



(A65) 



APPENDIX B: GAIN IN SPECTRAL DENSITY 
FOR SIMPLIFIED MODEL WITHOUT OPTICAL 
LOSSES 

Starting with this Appendix we simplify the model sys- 
tem. In particular, we assume here, that there is no op- 
tical losses in the system we examine. When considering 
the frequency-dependent rigidity based system, the ap- 
proximate formula (@J is used. 



Conventional harmonic oscillator 



2. Frequency-dependent rigidity 

a. Double-resonance case 

Consider now an oscillator with the frequency- 
dependent rigidity Q). In this case, 



2_ m^{Ql~nYi^--^^? 



(BIO) 



[see Eq. 0]. 

If the double-resonance condition © is fulfilled, then 
in close vicinity of the resonance frequency Qq = 



Smin"0 



(Bll) 
(B12) 



Performing again the same optimization as in previous 
subsection, we can obtain that again S,-min — ^o/\/2, and 



(B13) 



Spectral density of the total meter noise (|37|l is equal 



to: 



s,^^{^)=^SFm + \z{n)\'s^{n), (bi) 



where Z{yt) is the spectral image of the differential op- 
erator Z. In case of a conventional harmonic oscillator, 



(B2) 



and in close vicinity of the resonance frequency fio 



and 



where 



em 



SsuTSl (^) 



Smm ' 



= — iio , <c f^o , 

^2 _ Sf{^q) 



Sn 



2hmQ, 



2 ' 




(B3) 

(B4) 

(B5) 
(B6) 



see Eqsl P lI ^ . 

Let us require that £ (O) does not exceed a given value 
£o within as wide a frequency band AQ as possible. It is 
easy to show that this requirement is met if 



and 



Therefore, 



AO = 2£f„i„0o 



2 _ AO 



(B7) 
(B8) 
(B9) 



b. Two close resonances 



In the sub critical pumping case fiOfl . 

|Z(0)|2w 777,2(41/2 -Ogr;2)2, (B14) 



r(f^)«Cin + 



2 , (47.2-^^2^2)2 



(B15) 



These functions has a local maximum at 7/ = and two 
minima at v = ±Oo?y/2. 

The same optimization as in two previous cases gives, 
that 



(B16) 
(B17) 



£(0o) = Co - V2^nun , 

1 AO 

and the optimal value of parameter 77 is equal to 

Vc^Co, (B18) 



APPENDIX C: GAIN IN SIGNAL-TO-NOISE 
RATIO FOR SIMPLIFIED MODEL WITHOUT 
OPTICAL LOSSES 

1. Free test masses interferometer 

Rewrite the signal-to-noise ratio (|24|l as follows: 



SNR = - 



|j^signal(»)pdO 
'5'sum(^) 



(Cl) 



For conventional interferometer (without optical 
springs), the total meter noise spectral density is equal 
to (see 



20^ 



0^(02 + 02) 



02-^02 



20* 



(C2) 



15 



Therefore, in this case, 



3. Frequency-dependent rigidity 



SNR: 



\h{n)\^dn 



|fe(r>o)pr>o 



where 



\h[VL)\^dn 



, (C3) 



(C4) 



is the numeric factor depending on the gravitation-wave 
signal shape h{Vl) /\h{Q,Q)\. 



1. Conventional harmonic oscillator 

Substituting Eqs. (HH) and into Eq. (|rTT|l . we ob- 

tain, that for a conventional harmonic oscillator the 
signal-to-noise ratio is equal to: 



SNR = 



I Jsignal(^o)P 



TT Jo SF{n) + m^s^{n){ni-n'^y 



(C5) 



In the narrow-band case [see Eqs. H43LIB5P ]. this equation 
can be presented as the following: 



du 



^SOl(^o) 



(C6) 



In similar way, using Eqs. (|B1IB14() and (|C1|I . we ob- 
tain, that for optical spring based oscillator. 



SNR« -|F3ignal(r!o)|^ 
TT 



dv 



(C7) 



_ 72|fe(no)p»o Co 

In a pure double resonance case (ry = 0), 



Slightly better result can be obtained for the case of two 
optimally placed resonances. If 



^ ^/3 " 



(C9) 



then 



If the separation between the two resonance frequen- 
cies is too big, T] ^ ^0 (but still 77 ^ 1), then 



(Cll) 
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